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1.  Introduction. 


It  has  been  observed  over  years  that  experimental  relaxation  behaviour 

in  many  glassy  materials  is  generally  well  fitted  with  Wi lliams-Watts  function 

4>(t)  =  exp(-t/x  )a,  0  <  a  <  1  (1.1) 

c 

with  a  and  ie  constants  for  a  given  material,  where  a  is  an  index  of  slowness 

and  t  is  an  effective  relaxation  time.  This  is  in  contrast  to  the  conventional 
e 

Debye  exponential  form 

4>(t)  =  exp(-t/TQ) 

where  a  =  1  and  is  the  primitive  (Debye)  relaxation  time.  Such  stretched 
exponential  behaviour  has  been  noted  in  a  wide  variety  of  materials  including 
organic  liquids  and  solutions,  polymers,  and  in  particular  in  many  glassy 
materials.  The  recent  interest  in  relaxation  problem^5  ^  provides  a  number 
of  models  explaining  the  universality  of  formula  (1.1).  However,  a  connection 
with,  the  well  known  in  probability  theory,  class  of  Levy  a-stable  distribu¬ 
tions^1'^  is  not  realized,  except  the  very  recent  work  of  Montroll  and  Bendler11. 

Many  scaling  relations  for  complex  systems  in  the  physical  sciences  involve 
non-integer  exponents,  similarly  as  the  Fourier  transforms  of  Levy  a-stable 
distributions.  Consequently,  there  are  many  results  in  statistical  physics 

where  Levy  a-stable  distributions  are  invoked,  but  they  arc  not  recognized  or 

(14) 

kept  hidden.  For  example,  Khalfin  derived  an  asymptotic  formula  for  quantum 
description  of  decaying  physical  system  which  is  the  Fourier  transform  of  Levy 
a-stable  distribution  of  the  system  energy  with  0<a<l.  His  paper  influenced  many 
physical  studies  in  the  decay  theory  of  quantum  systems  but  this  relation  was 

not  recognized.  In  the  theory  of  dielectric  relaxation  of  polymeric  systems  and 

(3) 

glasses  Williams  et  al.  have  derived  a  series  representation  of  the  normalized 


2 


dielectric  parameter.  Its  imaginary  part  up  to  multiplication  factor  is  nothing 
else  as  known  series  representation  of  Levy  a-stable  density.  It  is  also  well 
known  that  this  series  representation  has  poor  convergence  properties  for 
large  arguments.  Neverless  this  fact  is  rediscovered  again  by  Lindsey  and 
Patterson^  who  were  forced  to  resort  to  special  multiprecision  computer 
programs.  Even  at  that  they  obtain  only  numerical  results  although  the  functional 

form  of  the  distribution  function  for  large  relaxation  times  is  of  equal  interest. 

f  71 

Helfand'-  continued  their  work  and  derived  an  asymptotic  formula  for  the  distri¬ 
bution  function  for  large  relaxation  times  which  gives  some  correction  terms  to 
the,  known  for  probabilists ,  series  expansion  of  Levy  a-stable  densities^^  .  He 
also  made  use  of  Zolotarev^lfl^  explicit  expressions  of  these  densities  for 
a  =  1/3  and  a  =  2/3. 

In  our  opinion  a  connection  between  the  statistical  aspects  of  Williams- 

Watts  relaxation  model  and  the  class  of  L£vy  a-stable  distributions  is  not  only 

a  coincident  analogy,  but  it  seems  to  capture  the  essential  statistical  nature 

of  relaxation  in  complex  materials.  Recall,  that  L6vy  a-stable  distributions 

form  a  class  of  universal  limit  laws  in  the  sense  that  they  describe  completely 

the  collective  i.e.,  macroscopic  behaviour  of  a  complex  system  expressed  by  any 

normalized  sum  of  independent  identically  distributed  random  quantities.  As  a 

consequence  of  this  important  theorem  in  probability^^’  ^  we  shall  derive  in 

this  paper  the  Williams-Watts  form  (1.1)  of  the  relaxation  function,  as  well  as 

the  existence  of  second  universality  relation,  first  obtained  by  Ngai^, 

connecting  i  with  t_. 

e  0 

Ngai's  cooperative-relaxation  model  based  on  Wigner's  random  matrices  theory 
of  energy- level  structures  leads  to  the  same  formula  (1.1)  and  the  effective 
relaxation  time  (in  his  notation)  is  related  to  the  primitive  relaxation 
time  i  as  follows: 

=  (l-n)(eY°w  )"  T  ,  0  <  n  <  1,  (1.3) 


where  y°  -  0.5772  is  Euler' s  constant  and  is  the  high-frequency  cutoff  of 

f  1 7") 

the  linear  density  of  excitations.  Recently,  Chamberlin  et  al.  have 
proposed  a  different  empirical  relation  in  the  context  of  the  remanent 
magnetization  in  spin-glasses 


1  —  n  . .  .  ...  1  —  n 

tp  =  (l-n)/L  to 


0  <  n  <  1 , 


(1.4) 


where  C  and  co  are  two  temperature- independent  constants  throughout  the  spin- 

(17  181 

glass  region.  The  controversy  which  arises  herev  *  ’  motivates  us  to  seek 

this  relation  from  purely  statistical  point  of  view  by  using  the  scaling 
property  of  the  characteristic  function  of  completely  asymmetric  Levy  a-stable 
distributions  relevant  to  the  relaxation  rates. 

Since  the  statistical  aspects  of  Williams-Watts  relaxation  model  are 
directly  connected  with  the  L6vy  a-stable  distributions  we  employ  them  to  study 
the  distribution  function  of  relaxation  times  which  complements  the  discussion 
given  by  Lindsey  and  Patterson and  Helfand^.  Also  the  statistical 
properties  of  dielectric  relaxation  are  discussed  from  this  same  point  of  view. 


|  Access' r~>  ‘"',T 
r  v-rr  l  ’  i  T 

|  ii  A  i  - 
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2.  Levy  ot-Stable  Distributions 

In  this  section  we  collect  some  basic  facts  on  Levy  ot-stabic  (LaS  in  short 

(13  19  20) 

notation)  distributions  needed  for  the  purpose  of  this  article  *  ’  .  One 

should  not  be  surprised  that  probabilistic  ideas,  in  particular  limit  theorems 
which  provide  a  macroscopic  distribution  for  a  collection  of  microscopic 
quantities,  play  a  role  in  statistical  physics.  We  can  refer,  for  example  to 
Khinchin' s^21^  treatment  of  statistical  mechanics  which  relics  heavily  on  so 
called  central  limit  theorem  leading  to  Gaussian  limit  distribution.  In  recent 

years,  inverse  power  long  tails  have  become  more  evident  in  the  analysis  of 

(22) 

physical  phenomena  of  condensed  matter  and  therefore  LaS  distributions  pro¬ 
vide  useful  models.  However,  the  principal  feature  of  these  distribution,  that 
we  want  to  stress  for  our  purpose,  is  that  they  arc  completely  described  as 
limits  of  the  normalized  sums  of  independent  and  identically  distributed 
summands.  Consequently,  LaS  distributions  express  some  kind  of  universality 
law. 

The  distribution  function  F(x)  is  called  stable  if  for  every  a^  >  0, 
bj ,  a^  >  0,  b2  there  correspond  constants  a  >  0 ,  b  such  that  the  equation 

FfajX  +  b  )  *  l(a2x  +  b2)  =  F(ax  ♦  b)  (2.1) 

holds.  The  symbol  F^  *  l'2  indicates  the  convolution  of  two  distributions  in 

the  sense 

Fj  *  l'2  =  J  FjU-y)  dl-2(y).  (2.2) 

It  turns  out  that  always 

a  =  (a^  +  a“)  1A*  with  0  <  a  <  2  (2.3) 

and  the  constant  a  is  called  the  characteristic  exponent  of  LaS-distribut ion. 


5 


Equation  (2.1)  can  be  solved  in  terms  of  characteristic  functions,  i.e.,  via 


Fourier  transform 


f(s)  =  exp[ isx]dF(x) . 


Thus  for  the  distribution  function  F(x)  to  be  LaS  it  is  necessary  and 
(19) 

sufficient  tliat  its  characteristic  function  f(s)  be  represented  by  the 
f ormul a 


iys  -  o  |  s  |  fl  -  iB  s  Lgn(s)tan(7ta/2) }  if  a/1 
log  f(s)  =  (2.4) 

iys  -  o|s|  -  iB(2/Ti)os  log  |s|  if  a  =  1 


where  a,B,y  and  a  are  real  constants  with  a  >  0,  0  <  u  <  2  and  |3|  <  1.  Here 
(t  is  the  characteristic  exponent,  y  and  o  determine  location  and  scale.  The 
coefficient  3  indicates  whether  the  LaS-distribution  is  symmetric  (B=0)  or 
completely  asymmetric  (]B|  *  1). 

Since  f(s)  is  absolutely  integrable  the  corresponding  LaS-distribution 
has  a  density.  It  will  be  denoted  by  p(x;a,B)  and  can  be  calculated  from  the 
Fourier  inversion  formula.  It  seems  impossible  to  express  them  in  a  closed 
form,  but  series  expansions  arc  available.  Only  in  the  case  0  <  a  <  1  the 
laS-densitics  with  |B|  =  1  are  one-sided  i.e.,  their  support  is  [o,+°°)  for 
B  =  1  and  (-m,()]  for  p  =  -1. 

By  means  of  a  linear  transformation  of  the  variable  in  the  distribution 
function  F(x),  it  is  possible  to  make  0=1  and  y  =  0.  Observe  that  only  the 
case  a  >  ()  is  of  interest,  because  o  =  0  corresponds  to  a  degenerate  distri¬ 
bution.  We  will  use  the  following  asymptotic  formulas^15)  for  LaS-densities 
when  0  <  a  <  1 ,  y  =  0  and  o  =  1 . 


P(x;ot, B)  =  (l/7i)  y  A  x 
n=l  n 


-na-1 


for  x  -» 


(2.5) 


.  (-1)  T(na+1)  ,,  a2 .  2  na.n/2  .  rtiu  tiu.t 

=  - — — pp — - -  (1+8  tan  j- )  sm  nLy-  +  a  ret  an  (8  tan— ^-)J 


and 


p(x;a,l)  = 


B(a)  A (a) 


2iT(l-a) 


-A(a) 


-1 


x  2  exp(-B(a)x"^a^ ) 


[l  +  0 


A  (a) 


-  e 


]  for  x  -+  0 , 


where 


and 


a 

B(a)  =  (l-g)a1_a  (cos  ™) 


1-a 


A  (a) 


a 


1-a 


(2.6) 

(2.7) 

(2.8) 


The  most  convenient  formulation  of  the  limit  theorem  which  gives  the 
description  of  the  distribution  law  F  governing  the  sum  of  a  large  number  of 
mutually  independent  identically  distributed  random  quantities  X^  can  be  given 
in  the  following  form^*^ :  Only  stable  distributions  have  a  domain  of 
attraction ,  i.e.,  there  exist  norming  constant  an  >  0,  b^  such  that  the  distri¬ 
bution  of  a_1(X.  +  X_  +...+  X  )-b  tends  to  F.  The  normalization  constants  can 
n  v  1  2  n  n 

be  chosen  in  such  a  way  that  a  =  n^a. 

'  n 

If  ('.  denotes  a  common  distribution  of  the  random  variables  X^  put  for 
x  >  0, 


u(x) 


y2C(dy). 

“X 


(2.9) 


Then  in  order  that  G  belongs  to  some  domain  of  attraction  it  is  necessary  that 

2  —  Q, 

u  be  of  the  form  u(x)  ~  x  S(x),  where  0  <  a  <  2  and  S  is  slowly  varying  at 
<»,  that  is,  for  fixed  x  >  0  S(tx)/S(t)  -+  1  when  t  •+  °».  There  arc  only  two 


cases : 


(ij  i  f  nt  =  2  then  G  belongs  to  the  domain  of  attraction  of  Gaussian 


distribution  if  and  only  if  u  is  slowly  varying, 

(ii)  if  0  <  n  <  2  then  G  belongs  to  the  domain  of  attraction  of  LaS- 

2 -ct 

distribution  if  and  only  if  u(x)  ~  x  S(x)  and  the  following  two 
limits  exist 


__  1  -  (:(x) 

1  -  G (x)  ♦  G ( x ) 


_ C(-x) 

1  -  G(x)  +  G(-x) 


(.2.10 


Observe  that  case  (i )  holds  whenever  the  second  moment  of  G  exists,  i.e., 
Var  G  *"  «>.  for  di st ri butions  without  variance  we  have  case  (ii). 


3.  The  Wi 1 1 iams-Watts  Function. 


One  of  the  most  characteristic  features  of  glassy  materials  is  the  non¬ 
exponential  character  of  relaxation  as  a  function  of  time.  This  has  been  shown 

(231 

to  be  true  in  a  wide  variety  of  materials,  including  organic  liquids  and 
r  '>4)  (4  5  251 

solutions  “  ,  polymers  ’  ’  ,  and  in  particular,  many  glass  forming  material 

It  is  a  striking  fact  that,  despite  the  variety  of  materials  used  and  of  the 
experimental  techniques  employed,  the  relaxation  behaviour  is  universal^. 

Good  fits  to  the  observations  have  been  obtained  with  a  functional  form  termed 
the  Wi 1 1 iams-Watts  function  given  by  formula  (1.1). 

It  is  common^’ ’ ^ ^  to  attempt  to  interpret  general  relaxation  in  terms 
of  exponential  relaxations,  writing 


-t/t  , 
o(t 


4>(t)  =  e  p( 
J0 


)  dx, 


where  p(x)  is  the  density  of  a  statistical  distribution,  i.c.,  p(x)  dT  =  1 , 

of  relaxation  times  x  across  different  atoms,  clusters,  or  degrees  of  freedom. 
If  u  =  x  /x  where  x^  is  a  single  relevant  relaxation  time  associated  with 
Debye  relaxation,  then  u  is  called  a  relaxation  rate  and  is  intcrprctable  as 


dimensionless  time.  Substituting  s  =  t/x^  in  formula  (3.1)  we  have 


■*  S  /  U  -  * 

<1>(X  s)  =  e  '  X  U  Xx (./u)du 

j0 


-sn 

-  e  w 
0 


(U)dp, 


(3.2) 


where 


1  (b)  =  T  u’  p(T_/u) 


(3.3) 


is  the  density  of  a  distribution  of  dimensionless  relaxation  rates.  Since  this 
approach  is  microscopically  arbitrary  we  may  consider  the  random  variables 


p  =  r  /i  .is  the  possible  relaxation  rates  of  elements  in  a  given  complex 
material.  The  index  i  indicates  the  number  of  the  elements  in  the  system 
and  p.  are  mutually  independent  and  identically  distributed.  The  problem  is 
then  to  determine  a  random  variable  p  for  describing  the  collective,  i.e., 
macroscopic,  behaviour  of  the  relaxing  system.  The  idea  is  to  introduce  a 
suitable  normalized  block-variable. 

Suppose  wc  divide  our  system  in  blocks  of  length  n.  The  total  relaxa- 

n 

tion  rate  of  each  block  will  be  equal  to  £  p. .  Clearly  this  is  not  a  good 

i  =  l  1 

collective  variable  since  for  large  n  it  goes  to  infinity.  Therefore  wc  need 
to  consider  normalization 


„  =  „i. 

i  =  1 


(3.4) 


Thus  it  follows  from  the  limit  theorem  discussed  in  Section  2  that  the  distri¬ 
bution  w(P)dp  of  the  collective,  i.e.,  macroscopic  relaxation  rate,  belongs 
to  the  class  of  LaS-d i stri but  ions  0«t<2.  Since  relaxation  rates  are  non-nega- 
tive,  w(p)  has  to  have  only  non-negative  support.  Consequently,  w(p)  is  the 
I,uS  density  with  ()<u<l  and  3=1,  (cf.  Section  2), 


w(p)  =  p(p;u,  1)  ,  ()<u<l. 


13.5) 


Wc  have  then  shown  that  the  macroscopic  behaviour  of  relaxation  rates  of 
our  system  is  described  by  a  distribution  belonging  to  a  rather  restricted 
class  of  l.rtS-d  i  stri  but  ions  with  parameters  0<cx<l  and  8=1.  The  possible  initial 
distributions  of  the  relaxation  rates  split  naturally  into  universality  classes 
parametrized  by  the  characteristic  exponent  u,  according  to  which  domain  of 
attraction  they  belong.  Observe  that  all  relaxation  rates  have  an  infinite 
variance  and  mean  value,  since  for  I/tS  distributions  only  p-th  moments  for 


Now  using  (3.2)  we  get 


4>(t0s) 


GO 

e~SU  p(ii;a,l)dM, 

0 


(3.6) 


i.e.,  4>(tqs)  is  the  Laplace  transform  of  LaS  density  p(p;a,l).  Recall,  that 
from  (2.4)  the  Fourier  transform  of  p(p;a,l)  can  he  written  as 


f(s)  =  exp  {iys  -  os  (1  -  i  tan(m/2) ) } 

(X 

=  exp  {iys  -  (o/cos  (tra/2)  )s  (cos(tra/2)  -  i  sin(TT«/2))} 

=  exp  {iys  -  [o/cos (tra/2) )  (-is)u)  .  (3.7) 


By  using  n  general  relation  between  Laplace  and  Fourier  transform  for  functions 
with  the  non-negative  support  L(p(y);s)  =  F(p(ji);is),  (3.6)  and  (3.7)  give 

<t>(x(1s)  =  exp  {-ys  -  (o/cos  (tra/2))  s®}  .  (3.8) 

Since  in  our  notation  is  =  t  thus  (3.8)  can  he  rewritten  as 

<j>(t)  =  exp  {-yt/T()  -  (o/cos  (na/2))  (t/x0)a}  (3.9) 

As  we  observed  in  Section  2  one  can  choose  y  =  0  and  hence 

<)>(t)  =  exp  f-(t/Tc)a},  (3.10) 

where 

t®  =  f  ( I/O)  cos  (ticx/2)  |  Tq  .  (3.11) 

Thus  (3.10)  gives  the  Wi 1 1 iams-Watts  form  of  the  relaxation  function  and  (3.11) 

gives  the  relation  between  the  effective  relaxation  time  Tc  and  the  primitive 

relaxation  time  t^.  Let's  remark  that  formulas  (3.10)  and  (3.11)  have  universal 

character  as  a  consequence  of  the  use  of  universal  limit  law  in  macroscopic 

behaviour  of  the  relaxing  system.  1111 s  provides  an  answer  to  the  question 

(12) 

which  has  been  posed  recently  by  Palmer  ct  al 


Since  cos  (Tia/2)  ~  1  -  n  for  ()<ot<l  then  (3.11)  can  be  approximated  by 


4.  The  distribution  Function  of  Relaxation  Times 


If  the  statistical  approach  to  explaining  nonexponential  relaxation 
behaviour  is  taken,  it  becomes  important  to  determine  the  distribution  of 
exponential  relaxation  times  which  produce  the  Wi  1 1  iams-Watts  function^'7, *  ^  _ 
In  this  section  we  will  employ  the  known  asymptotic  formulas  for  one-sided 
LaS-dens i t ies  to  derive  the  asymptotic  behaviour  of  the  distribution  function 

(in  probabilistic  terminology  density)  p(t)  ,  defined  in  (3.1). 

- 1  7 

Since  p(t)  =  t  p  w(p)  according  to  (3.3),  thus  by  (2.5),  if  we  put  3=1, 
we  have  for  i  -*■  0 


0(T)  = 


(-1)  T  (not 


1) 


TIT 


e  n=l 


n ! 


an  - 1 


stn  nira. 


(4.1) 


where 


Tc  =  |cos(TTa/2)]1/aT() 


For  t  -*■  "•  we  get  from  (2.0) 

A  (a) 


P(t)  = 


T  /Imx 
c 


4  >(«)-! 


1 

• 

’  ' 

A  (a) 

T 

exp  ■ 

" 

T 

c 

(4.2) 


where  i  =  B(a)~^  and  arc  8ivcn  by  (2.7)  *',nd  (2.8), 

respect i vely. 

Observe  that  the  existence  of  average  relaxation  time  <t>,  and  higher 
moments,  <tn>,  n>l  follow  without  any  trouble  since  by  (3.3)  and  (2.5) 


n 

<T  >  = 


n  ,  .  ,  n 
T  0(T)dt  =  Tq 


U~nw(p)dp 


(4.3) 


-a-l-n 

P  db  <  00 


This  is  in  sharp  contrast  to  the  moment  less  distribution  of  the  relaxation 
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rates  ii  =  i  /t. 

The  relationship  between  the  relaxation  function  <J>{t)  and  the  LaS- 
ilensitv  w(p),  sec  (3.2),  permits  to  find  the  relaxation  time  <t>  and  higher  moments 
without  explicit  knowledge  of  p(x).  Namely,  since 


w  (p) 


2  iii 


i°° 


-joo 


e,JS  4>('t()s)  ds. 


where  s  =  t/i()  and  p  =  one  has  from  (4.3) 


<rn>  «  rjj 


,  -n  ,  1 

U  ^  2rri 


i'ion 


-i°° 


^  Ct0s) ds 1  dp. 


Moving  the  p  ”  term  inside  the  integral  over  ds ,  interchanging  the  order  of 
integration,  and  moving  the  <|>(t0s)  term  outside  of  the  integration  over  dp 
yields 


<-  rn>  s  TJj 


♦(T0#)  l2Tii 


1  f10°  -n  pS 

r  J  p  e 


dp]  ds 


■i°° 


{ T^/I’Cn)  > 


<f»(T()s)  Sn_1ds  =  (l/r(n)} 


4>(t)  tn_  1  dt 


=  f l/r(n) } 


..n-1  t  I  l\a\  a* 
t  exp  l - (— — )  )  dt 

0  e 


t"  r(5)/(or(n)) 


t"  [(I/O)  cos  (TO/2)]n/0t  r(J)/(«r(n)), 


(4.4) 

(4.5) 


where  we  have  used  formulas  (3.10)  and  (3.11).  formula  (4.4)  corresponds  to 
the  result  in  l.indsey  and  Patterson ^  . 


.  ‘  j 

.  *  .  L.  *  •  -  .  , 


-'.V.-V/.-.V/,.’.'. 


’.v  /. 


5.  Dielectric  Relaxation 


The  complex  dielectric  permittivity  e* (u>)  is  related  to  the  relaxation 
function  <J)(t)  by  a  one-sided  Fourier  transform^ ^ 


e*(w) 


0 


-  e 


/OO 

-iiiit 

f  d 

e 

0 

['dt  *lt) 

dt, 


(5.1) 


where  l-  (d/ dt) <p C t ) |  is  the  normalized  transient  decay  current  which  follows 
when  a  steady  electric  field  is  removed  from  a  sample.  Here  is  the  static 
dielectrical  constant  and  the  high-frequency  limit  of  the  dielectric 
constant.  One  generally  writes 


g  * (w)  -e^ 

— — — °°  =  e'(w)  -  ie"(w)  ,  (5.2) 

0  <*> 


where  e'(w)  and  f."(u))  are,  respectively,  the  real  and  imaginary  parts  of  the 
normalized  dielectric  parameter.  Williams  et  al.  derived  an  expression  for 
e*(w)  when  4>(t)  =  exp  {-(t/Te)a}  in  the  form 

g’(w)  -  ie"(u>)  =  l  (-l)n_1  - - —  r(-na  |  11  x 

n=l  (un  )"“  ”• 

e  (5.5 

,  tt2  .  .  tt2. 

x  (cos  n  ^ - 1  sin  n  j-) . 

By  the  comparison  with  (2.5)  one  obtains  immediately  that 

g" (a))  =  xiTp(x;a,())  ,  (5.4) 

where  x  =  wtc.  This  direct  relation  between  the  imaginary  part  r"  of  the 
normalized  dielectric  parameter  and  the  symmetric  1/xS  density  p(x;n,0)  has 


It  is  interesting  to 


been  ;ilso  shown  recently  by  a  different  method1  ; 
explain  why  here  in  contrast  to  Sections  3  and  4  symmetric  (B=0)  I,aS  distri¬ 
butions  arise. 

lor  this  observe  that  i f  we  put  y  =  0  and  o  =  1  in  formula  (2.4)  then 
for  0  <  (x  <  1 

log  f (s')  =  -  |s|a  (l  -  iB  sign  (s)  tan  (-na/2)}.  (5.5) 

Knowing  that  the  density  p(x;a,B)  of  such  LaS  distribution  is  real,  it  can  be 
written,  after  a  short  calculation,  in  the  form 

p(x;u,3)  =  Re  e  lsxf(s)  ds  = 

_oo 

sa 

c  cos  fxs-Bsa  tan(ira/2)]  ds.  (5.6) 

On  the  other  hand  the  integration  by  parts  in  (5.1)  yields 

-oo  a  j-°o  a 

[e'(u))  -  1|  -  ie"(ui)  =  -x  e  S  sin(sx)dx  -  ix  e  s  cos(sx)ds,  (5.7) 

■*0  ^0 

where  s  =  t/r^  and  x  =  mt  .  The  comparison  of  (5.6)  and  (5.7)  shows  that  the 
imaginary  part  t "(w)  of  the  normalized  dielectric  parameter  is  related  to 
LaS  distribution  p(x;a,B)  if  and  only  if  B  =  0,  i.c.,  when  the  LaS  distribution  is 
a  symmetric  one.  Namely  one  obtains  from  (2.5)  and  (3.11) 

f"(u>)  =  TiuiT^  p  (,,,Te  ;<*»()) 

v  ,  ,  .n+1  l’(na  +  1)  ,  . -na  .  ,  ,.,.,-n 

=  }.  (-U  - — i -  (wt  )  (cos(na/2)]  sin  (rnra/2) . 


(5.8) 


(>.  Concluding  Remarks 


There  are  several  approaches  to  explaining  nonexponential  relaxation 

( i_  i  ?) 

behaviour  in  glassy  materials  .  The  statistical  approach  interprets 

the  nonexponential  relaxation  behaviour  of  the  material  in  terms  of  a  super¬ 
position  of  exponentially  relaxing  processes  which  then  leads  to  a  distribu¬ 
tion  of  relaxation  times.  In  this  paper  we  derived  the  Wi 1 1 iams -Watts  formula 
directly  from  the  limit  theorem  in  probability  theory  that  explains  the  univer¬ 
sality  of  this  formula,  as  well  as  provides  a  relation  between  the  effective 

relaxation  time  t  and  the  primitive  relaxation  time  t,  , .  This  procedure  is 
e  O' 

different  from  those  in  the  existing  literature  and  provides  the  main  departure 

from  all  the  previous  work  on  the  statistical  approach  to  the  relaxation  problem 

The  most  important  consequence  is  that  the  theory  of  Levy  m-stable  distributions 

can  be  directly  employed  to  the  study  of  the  problem. 

We  do  not  propose  that  the  statistical  approach  is  the  physically  more 

correct  one;  instead  we  discuss  the  mathematical  foundation  and  consequences 

of  this  approach  which  must  be  understood  when  physical  interpretation  is 

attached  to  the  distribution  functions  underlying  the  relaxation.  We  would 

f  28) 

like  to  conclude  the  paper  by  quotation  from  l:.P.  Wigner1  "The  first  point 

is  that  mathematical  concepts  turn  up  in  entirely  unexpected  connections. 
Moreover,  they  often  permit  an  unexpectedly  close  and  accurate  description  of 
the  phenomena  in  these  connections.  Secondly,  just  because  of  this  circum¬ 
stance,  and  because  we  do  not  understand  the  reasons  of  their  usefulness,  we 
cannot  know  whether  a  theory  formulated  in  terms  of  mathematical  concepts  is 
uniquely  appropriate."  Let  it  be  an  excuse  for  the  authors,  if  they  were  not 
sufficiently  modest  in  some  formulation  of  this  paper. 
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